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Abstract 



In this note, by constructing suitable approximate solutions, we prove the existence of 
global weak solutions to the compressible Navier-Stokes equations with density-dependent 
' viscosity coefficients in the whole space M. N , N > 2 (or exterior domain), when the initial 

CNI . data are spherically symmetric. In particular, we prove the existence of spherically symmetric 

solutions to the Saint- Venant model for shallow water in the whole space (or exterior domain) . 
p j Keywords: Compressible Navier-Stokes equations; density-dependent viscosity coefficients 

< 

■ 1 Introduction 

+-> 

In this note, we consider the following compressible Navier-Stokes equations with density- 
dependent viscosity coefficients 

p t + div(pU) = 0, (1.1) 

> 
CO 

where t G (0, +oo) and x G R N , N > 2, p(x, t), U(x, t) and P(p) = p"> (7 > 1) stand for the fluid 
density, velocity and pressure respectively, 



(pU)t + div(pU ® U) - div(2%)£>(U)) - V(s(p)divU) + VP(p) = 0, (1.2) 



D(U) = i(VU + (VU) T ) 
is the strain tensor, h(p) and g(p) are the Lame viscosity coefficients satisfying 



>< ' h(p) > 0, 2h(p) + Ng(p) > 0. (1.3) 



In the last several decades, significant progress on the system (|l.ip - (|1.2p with positive con- 
stant viscosity coefficients has been achieved by many authors. In the case that the initial data 
are sufficiently regular and the initial density is bounded away from zero, there exists a unique 
local strong solution, and the solution exists globally in time provided that the initial data are 
small perturbations of an uniform non-vacuum state. For details, we refer the readers to papers 
[20] and the references therein. The situation becomes more complex in the general case of 
nonnegative initial density, and a number of important questions are still open. For example, 
the uniqueness of global weak solutions. The first general result on the existence of global weak 
solutions was obtained by Lions in [18]. There have been many generalizations of this result, 
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see [HI E31 [151 HU [23] . Using the compatibility condition, Salvi-Straskraba [22] and Choe-Kim 
[6] obtained the existence and uniqueness of the local strong solution. 

The results in [1 1 1. \19\ EI] show that the compressible Navier-Stokes system with constant 
viscosity coefficients have the singularity in the presence of vacuum. By some physical consid- 
erations, Liu, Xin and Yang in [19] introduced the modified Navier-Stokes system with density- 
dependent viscosity coefficients. As remarked in [19] . in the derivation of the Navier-Stokes 
equations from the Boltzmann equation through the Chapman-Enskog expansion to the second 
order, the viscosity is a function of the temperature, and correspondingly depends on the density 
for isentropic fluids. Meanwhile, in geophysical flows, many mathematical models correspond 
to (jl.ip - (jl.2p . In particular, the viscous Saint- Venant system for shallow water is expressed 
exactly as ([H|) - (Q]I with N = 2, h(p) = p, g(p) = and 7 = 2 (pQ El US]). As remarked 
in [10], new mathematical challenges are encountered for the shallow water equations and the 
multi-dimensional compressible Navier-Stokes equations (|l.l|) - (|1.2p . The main difficulty is that 
the velocity can not be defined in the vacuum state. 

For one-dimensional compressible Navier-Stokes equations (|l.ip - (|1.2p with h(p) = p e and 
g{p) = 0, 9 £ (0,1), there are many literatures on the well-posedness theory of the solutions, 
see [U [T^l [T71 [19l [2"5| ES] • Considering the free boundary problem of the spherically symmetric 
system, the local existence and uniqueness of the weak solution were obtain in [5], the large-time 
behavior of the global solution for data close to equilibrium was obtained in |27[ I28j . However, 
few results are available for multi-dimensional problems. In [JJ, Bresch, Desjardins and Lin 
showed the existence of global weak solutions in dimension 2 or 3 for the Korteweg's system 
with the Korteweg stress tensor kpVAp. An interesting new entropy estimate is established in [T] 
in a priori way, which provided some high regularity for the density. Later, a similar result was 
obtained in [2] with an additional quadratic friction term rp|U|U. Recently, Mellet and Vasseur 
[2T] proved the L 1 stability of weak solutions of the system (|l.ip - (|1.2p with N = 2, 3 and 7 > 1, 
based on the new entropy estimate, extending the results in [HE] to the case r = k = 0. Bresch 
and Desjardins constructed approximate solutions for the viscous shallow water system with 
drag terms or capillarity term and for the compressible Navier-Stokes equations with the cold 
pressure in [3] , and proved the global existence of weak solutions to these systems in [31 0] . In 
|10j . Guo, Jiu and Xin constructed a class of approximate solutions and proved the existence of 
global weak solutions for the spherically symmetric compressible Navier-Stokes equations with 
density-dependent viscosity in a bounded domain (iV = 2, 3, 7 > 1). 

In this note, we will construct a class of approximate solutions and prove the global ex- 
istence of weak solutions for the spherically symmetric compressible Navier-Stokes equations 
with density-dependent viscosity in the whole space or exterior domain (N > 2, 7 > 1). Us- 
ing the method in [10], we can construct the approximate solutions on the annular domain 
{e < |x| < R} by solving the approximate systems of (|l.ip - (|1.2p with h e {p) = h(p) + ep e and 
g e (p) = g(p) + (0 — l)ep e instead of h(p) and g(p). Then, using the usual zero extensions as 
in [12\ I13j. we can construct the approximate solutions on the entire domain W N . But, the en- 
tropy estimates of approximate solutions do not hold on the entire domain M> N , only hold on the 
annular domain. Using some techniques in Proposition 13.31 we can prove that Vy/p belongs to 
L°°(0, T; L 2 (W N )), so that the nonlinear diffusion terms in the definition of weak solutions will 
make sense. The extension method in [10], can preserves the uniform L°°(0, T; H 1 ^)) estimate 
of VP^j but seems not applicable to build approximate solutions in the whole space or exterior 
domain. 
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2 Statement of the results. 

The Cauchy problem of the compressible Navier-Stokes equations can be written as 

pt + div(pU) = 0, (2.1) 

(pU) t + div(pU <g> U) - div(2%)£>(U)) - V( 5 (p)divU) + VP(p) = 0, (2.2) 
with initial conditions 

p\t=o = po > 0, p\J\ t =o = mo- (2.3) 

Before introducing the notion of weak solution, let us state the assumptions on the viscosity 
coefficients, as in |21|. 



Conditions on h(p) and g(p): 

We assume that h(p) and g(p) are two C 2 (0,oo) functions satisfying 

g(p) = 2ph'(p)-2h(p), (2.4) 

h'{p) > u, h(0) > 0, (2.5) 

\ 9 '(p)\ < h'ip), (2.6) 

uih(p) < 2h(p) + Ng(p) < u 2 h(p), (2.7) 
where v G (0, 1) and zaj > ^i > are three constants satisfying 



AN - 4V2iV 2 - 4N + 4 ^ - 2 4iV + 4V2iV 2 - 4iV + 4 i/ 2 - 2 

iV 2 - 4iV + 4 < iV ' iV 2 — AN + 4 > N ' ~ ' ' 

When iV > 3 and 7 > t^,, we also require that 



lim inf w=r~ — > 0, 



for some small e > 0. 

Remark 2.1. From the above conditions, one has 



jV— 1 , v\ JV— 1 1 ^2 

Cp~ h 2i7 < /l(p) < Cp~ r 2JV) p > 1, 
N — l 1 ^2 JV— 1 1 ^1 

Cp" h 2F < < Cp~ h 2iv , p < 1. 



(2.10) 



Definition 2.1. We say that (p,U) is a weak solution of (pTT| - (pl3"]) on R N x [0,T], provided 
that 

(1) 

p G L^TjL 1 ni 7 (l^)), L^^T-H 1 ^)), 
JpUeL°°(0,T; (L 2 ^))^), 
%)D(U)GL 2 (0,T;(^ 1 ' 1 (R Ar )) AfxAr ), 5 (p)divU G L 2 (0, T; ^(M^)), 
with p > 0; 

(2) For any i 2 > h > and 0i G C^(R N x [0, 00)), the mass equation (|2.ip holds in the 
following sense: 

I pfadxft = / 2 / (p^i + pU • V<h)dxdt; (2.11) 
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(3) The following equality holds for all smooth test function <p2{t,x) £ (C 2 (M. N x [0,00)))^ 
with 4> 2 {T,-) = 0: 



f m • (/> 2 (0, x)dx + f [ (Vp(VpU) • ^<fe + VpU ® ^/pU : V</> 2 ) dxcft 

+ / / p^divfcdxdt- < 2/i(p)D(U), V0 2 > - < <?(p)divU, div(/> 2 >= 0, (2.12) 
Jo JR N 

where the diffusion terms make sense when written as 

<2%)£>(U),V0> = -/ ^(^pU^diifydxdt- [ (^pU j )2h'(p)d i ^pd i (l )j dxdt 

Jm. N \1P Jr n 

^/pUi)dji4>jdxdt - / (y/p[Ji)2ti '(p)djy/pdi</>jdx.dt, 



^fp 
Hp) i 

VP 

and 

< 5 (p)divU,div0 >= - / ^(y/pU^dijfadxdt - [ {Vf^ j )2g\ P )d j Vpdi(t>idxdt. 

Jr n VP 

In this paper, we will construct global spherically symmetric weak solutions to (I2,ip - (l2,3p , 
The initial data are assumed to satisfy 

po > a.e. in R N , m = a.e. on {x e R N \p (x) = 0}, (2.13) 

poetfntfiR"), ^6iV), ^(l + ln(l + ^|))eL 1 (IR iV ). (2.14) 
y/Po Po Po 



The main result of this paper is the following: 

Theorem 2.1. Assume that 7 > 1, h{p) and g(p) satisfy conditions {2.J$ - $2~9\i . If the initial 
data have the form 

Po = Po(|x|), m = m (|x|)- 

r 

and satisfy \2.13\) - §2.14% then the initial-value problem \2. ij) -( fiOp has a global spherically sym- 
metric weak solution 

p = p(I x M)> u = u(\x\,t)- 

r 

satisfying for all T > 0, 

KMjeCdO,! 1 ];! 1 ^))), (2.15) 



/ p(x,t)dx = / po(x)dx. (2.16) 
Jr n Jr n 



Moreover, it holds that 



sup / fi|V%)| 2 + P(i + |U| 2 )(l + ln(l + |U| 2 ))^x< c, (2.17) 
te[o,T]JR N \P J 

where C is a constant. 

Remark 2.2. Using the similar argument as that in [10], one can obtain that 

sup f p\V\ 2+v dx < C, 
te[o,T] Jr n 

when J RJV po\\Jo\ 2+v dx < C for some small 77 £ (0, 1). 



Remark 2.3. Similarly, using the usual zero extension method, one can obtain the similar result 
for the existence of global weak solutions for the spherically symmetric compressible Navier- 
Stokes equations with density-dependent viscosity in a bounded domain (N > 2, 7 > 1). 

Remark 2.4. Under conditions (|2.4p - (|2.9p . using the similar argument as that in [2.1], one can 
easily obtain the similar result as that in [21] with N > 2 and 7 > 1. 

• Exterior problem 

Using the similar proof of Theorem 12. 1[ we can study the following exterior problem: 

pt + div(pU) = 0, f>0,xe!], (2.18) 

{p\J) t + div(/>U ® U) - dw(2h(p)D{\J)) - V(s(p)divU) + VP(p) = 0, (2.19) 
with boundary and initial conditions 

(pU)| xean = 0, p\ t =o = Po>0, pV\ t =o = m , (2.20) 
where ll = {x£ R N \\x\ > 1}, N > 2. 



Definition 2.2. We say that (p,U) is a weak solution of ([2TT8|> - ([2T20]) on O x [0,T], provided 

(a) The condition (1) in Definition 12.11 where W N is replaced by fi; 

(b) For any i 2 > *i > and <pi £ C^(R N x [0, 00)), the mass equation (|2.ip holds in the 
following sense: 

f pfadxlll = [ 2 [ {pdt^i + pU • Vfa)dxdt; (2.21) 

JO, Jt! JU 

(c) The condition (3) in Definition 12.11 where M> N is replaced by f2. 



Using the similar proof of Theorem 12.11 and ||v / P( r )lli° ([i,oo)) ^ II vpWIIh 1 ([1,c»))j we can 
obtain the similar result without the condition (12. 9p . Here, we give the following theorem and 
omit the proof. 

Theorem 2.2. Assume that 7 > 1, h{p) and g(p) satisfy conditions \2.J$ - ^2~7fy . If the initial 
data have the form 

Po = Po(|x|), mo = m (|x|)- 

r 

and satisfy i2.13\) - (2.14\ ) where M N is replaced by Q, then the initial-value problem $2. l\) - ^2. &\) 
has a global spherically symmetric weak solution 

9 = P(I X M)> U = u(|x|,t)- 

r 

satisfying ^2~15\) - [2~17\ ) where R N is replaced by fl, for all T > 0. 

Remark 2.5. In particular, we get the existence of spherically symmetric solutions to the Saint- 
Venant model for shallow water system in the whole space or exterior domain. 
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3 Proof of Theorem [2J] 

The key point of the proof of Theorem l2,ll is to construct smooth approximate solutions satisfying 
the a priori estimates required in the L 1 stability analysis. The crucial issue is to obtain lower 
and upper bounds of the density. To this end, we study the following system as an approximate 
system of (l2TT]) -([2T2l. 

Pt + div(pU) = 0, (3.1) 

(pU) t + div(pU <g> U) - div((2%) + ep e )D{\])) 

-V((g(p) + {9- l)e/)divU) + VP(p) = 0, (3.2) 

where e > is a constant and 6 = N ~^ +a with a G (0, 1) satisfying 

T T , N , ,i/i-2 a-L ir , JV , v 2 -2 , , 

Vi(_) < m n{_, _}, V5(— ) > — , (3.3) 

where 



, , AN(m - 1) - 4JN 2 (m - l) 2 + (N - l)(m - l)(m - 2) 2 
Fl(m) = (iV_i) (m _2) 2 



and 



, , AN(m - 1) + 4A/iV 2 (m - l) 2 + (N - l)(m - l)(m - 2) 2 
V2(m) = 



(iV-l)(m-2) 2 

Remark 3.1. From (|2.8p . we can choose a small constant a satisfying ()3.3p . 
When p(x,i) = p(r,t), U(x,i) = u(r,t)f , (pn) - (pT2~]) becomes 

Pt + pu ) P + ^ ^- = 0, 3.4 

r 

+ puu r + ijP) r + (2/i + ep e ) ^ N ~ l)u = ((2h + g + 0ep e )(u r + (jV ~ 1) " )) r , (3.5) 

r r 

for r > 0. We will first construct the smooth solution of (|3.4p - (|3.5p in the truncated region 
0<£<r<i?<oo with the following boundary conditions and initial condition 

u(r, t)\ r=e = u(r, t)\ r=R = 0, (3.6) 

(p, u){r, 0) = {p ,e,R,5, u , s ,R,5) := (po, £ ,R * Js, ua,e,R *J S ), e < r < R, (3.7) 
where Js is a standard mollifier, 



Po,e,n( r ) 

and 

U0,£,R{ r ) = 



p (e) +£, r G [0,e], 
p (r) + e, r G [e,R], 
k p (R) +£, r G [JS,oo) 



0, rG[0,£ + 2<5], 

-^sg-, r€[e + 25,R-2S\, 
0, r e[R— 26, oo). 



We assume that e and -R satisfy eii^ < y^e. Letting e — ► and R —* oo, we can easily obtain 
that (po, £ ,r, u o,£,r) convergence to (po,uo) in spaces given in (|2.14|) . From (|3.3|) and similar 
arguments as that in [lOj, one can obtain the smooth solutions (p s ' R ' (r,t),u £ ' R ' (r,t)) to the 
approximate system (I3.4p - (l3.7p . 
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Remark 3.2. To obtain the existence of (p £ ' R ' S (r, t),u e ' R ' S (r, t)), we need to consider the following 
system in the Lagrangian coordinates: 

Pr + p 2 ^- 1 ^, = 0, 

r l ~ N u T + {pi) x = [(ph + pg + e0p e+1 ){r N - l u) x ] x -(h + ep e ) x {N - l> , 

r 

u(0,t) = u(1,t) = 0, 

(p,u)(-,0) = (po, £ ,R,6,UO,e,R,s)- 

N 

From (|3.3p and similar arguments as that in [10] . one can obtain that u G L°°(0, T; L x ~ a ), 

JV 

(/O e )x £ L°°(0,T; L x ~ a ) and G L^. (for simplicity, we omit the superscript). To estimate 
||ii|| jv , we need to estimate the following terms 

-(m - \){2hp + gp + e6p e+1 )r 2N - 2 u m ~ 2 ul 

-[2{N - l)ph + (N - l) 2 pg + e(N - 1)(0(N - 1) - N + 2)p 1+9 ]r- 2 p~ 2 u m 
-[pgm(N - 1) + em{N - 1)(0 - l)p 1+9 }r N ~ 2 p^u™' 1 ^, (3.8) 

where m = j^— . From (|3.3p . we have 

dSSD < -C(ph + ep e+l ){r 2N - 2 u m - 2 u 2 x + r^VV"). 

JV 

Then, using similar arguments as that in [10], one can obtain u G L°°(0, T; L x ~ a ). 
So far, 

are defined on e < r < i?. To take the limit (ej,Rj,8j) -> (0,oo,0), we extend (p 6 '^ 5 , w^'" 5 ) to 
the whole space in the following way 

^(r,.)-{f' A (3.9) 

s*m={ : ls 6 e . [e ' fl1 ' (3.10) 

For simplicity, we denote the obtained approximate solutions {p £ j' R j^3 ', u e i' Ri,Si } by n 3 '}. 
Let p>(x,t) = p j (r,t), \J j (x,t) = vP(r,t)f, B E>R = {x G M. N \ e < |x| < R} and B R = {x G 
|x| < i?}. 

Using similar arguments as that in proofs of Lemmas 3.2 and 4.1 in [10], and the similar 
argument as that in [18J (§5.5), we have the following lemma. 

Lemma 3.1. There exists a constant C independent of e, R and 5 such that 



sup / p>(x,t)dx<C, (3.11) 
te[o,T] Jr n 



1 1 



sup / ( ±pj\\]i\ 2 + — L — (pi)A (x,t)cfc + f T f (uxh(pi)\VW\ 2 )( X ,t)dxdt 
te[o,T]Ju N V 2 7~! / Jo Jr n 
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+ae [ [ ({p>) e \VXJ j \ 2 )(x,t)dxdt < C, if 7 > 1, (3.12) 
Jo Jr n 



sup / ( \p>\W\ 2 + pjlogfP -plog/9- (logp + l)(^ - p) J (x,i)dx 
te[o,T] Jk^ V z / 

+ / T / (n^')|VU 3 '| 2 )(x,t)cixcft + ae / T / ((^') e |VU J '| 2 )(x, t)dxdt 
Jo ii^ Jo Jr n 



'0 JlR^ JO JR^ 

< C, if 7 = 1, p = e~W, (3.13) 



/ 1 

JB e K 2 

e J '-"3 



2/>V) + tefpi^- 1 

e J '-"3 

1 / 2ti{p>) + ee(pi) 9 -\ 



sup / - p 1 

te[o,T] 



x, t)dx 



+ y y ( — vr y vr 7 — vp>v(p>y ) (x,t)d*dt < c. (3.14) 



Moreover, the following uniform estimate hold 

sup \\v^\\m(B E . R .) <C- (3-15) 

te[o,T] J ' J 

From this lemma, we can obtain the following lemma. 

Lemma 3.2. The pressure (/>?) 7 is bounded in (R N x [0, T]) when N > 3, in L@(R N x [0, T]) 
for all (3 G [1, 2) when N = 2. 

Proof. From (|23j> . ([3TTT ]) -(|3TH 1) . we have (/>0 7/2 is bounded in L 2 (0, T; H x {B £j:Rj )). 

When iV > 3, we get (z^) 7 / 2 is bounded in L 2 (0,T; (B £j:Rj )) or (/>?) 7 is bounded 

in L 1 (0,T;L^( J B £j;Rj )). Since (p?) 7 is bounded in L°°(0, T; L l {B EjAo )), Holder's inequality 

implies that (p^) 7 is bounded in L^h(B £>R . x [0,T]). From (|3.9p . we obtain that (p 3 ) 7 is 

bounded in L^(R N x [0,T]). 

Similarly, we can get that (p 7 ) 7 is bounded in L^(M. N x [0, T]) for all G [1)2) when 
N = 2. □ 

In the following proposition, we will estimate ||p|U| 2 ln(l + |U| 2 )||£i( R jv). 

Proposition 3.1. If v\h <2h + Ng < u 2 h and J °° p Q (l + |u | 2 ) ln(l + luo) 2 )^ 1 ^ < C, then 
the following estimate is true 



fRj \ u j\2 

/ p jl -^ln(l + \u : >\ 2 )r N - 1 dr <C (3.16) 
J £j 2 



sup 

*S[0,T] Je 



where C is a constant independent of 6j , Rj and Sj . 

Proof. Multiplying (|3.5|) by r N ~ x vP{l + ln(l + In-?! 2 )), integrating the resulting equation and 
using (|3.4p yield 



dt J £ . 2 
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+ I (2/ l + e( / >'y)(l + ln(l + |^| 2 ))((^) r 2 + (iV-l)^)r Ar - 1 dr 



R 



+ {2h + e{p>)») 



l + M 2 



R i 



+ (g+(6- l)e{p>Y)(l + ln(l + \v?\*)){{v?) r + (N- 1)— )V v - i (ir 



A', 



2(w 



+ / ((/^') 7 ) r (l + ln(l + H 2 )V r^dr = 0. 



Since v\h <2h + Ng < u 2 h and (1 + N(9 - l))e(p>) e = ae(f>>) e , we have 

d f R * 4 l + \u j \ 2 
dt 



ln(l + H 2 )^- 1 *- 



+ / (yxh + ae{p>Y){l + ln(l + \u j \ 2 ))((v?fi + (N - l)^J-) r »- l dr 



({p j V)r{l + ln(l + \v?\ 2 ))v?r N - l dr 



Ri 



< C I (h + e(p>) )((ui) 2 r + (N — l)^_)r N - l dr. 



(3.17) 



Using integration by parts and Young's inequality, we have 

r*R 

\{p>y) r {\ + ln(l + \u j \ 2 ))u j r N - l dr 



R i 



< C |^|(l + ln(l + |n J '| 2 ))(p J ') 7 r 7V - 1 dr + C / \u j \(l + \n(l + \u j \ 2 ))(p>yr N ~ 2 dr 



< 



i/] 1 ^ h(l + ln(l + H 2 ))((^') 2 + (N- l)^L)r N ~ l dr 



< 



-C I h- x { f P)^{l + \D.{l + \u i \ 2 ))r N - x dr 
" l f ' h(l + ln(l + \u J \ 2 )){(u J ) 2 + (N- \) { ^)r N - l dr 



+C 



R» 



l h 



£ \ 2-5 
2 \ 



2-5 
2 



r N-l dr 



R i 



^'(^^(l + l^'l 2 ))^^- 1 ^ . 



Combining it with (l3TT]) - (l3TT3| ) and d3T?]) . we get 

R < 1 + M 2 



sup / fP 

i£[0,T] Jej 



< C + C s 



Mi + i^'i 2 )^- 1 ^ 



-2. \ 2-5 



2-5 
2 



r N-l dr 



(3.18) 



From (|2.5p . we have h>vp and 
' ( (P 1 ) 



Then, using Lemma [3.21 we check that the right hand side is bounded L 1 in time for some small 
5, without any condition when N = 2, and when N > 3 under the condition that 

iV + 2 

which gives rise to the restriction 7 < J_ 2 - In either cases, we have 

rRi 1 + M 2 



/ + 1 in(i + i^'i 2 )^- 1 ^ < a 



sup 

te[0,T] 

When iV > 3 and 7 > , we need the extra hypothesis (|2.9p to show that the right hand 
side of (|3.18p is bounded and to obtain the same result. □ 

From ([3T9P - ([3TT0ll . we deduce that 
Corollary 3.1. 

r ITP'I 2 

sup / p> l —^ln(l + \\Ji\ 2 )dx<C. (3.19) 

JV AT 

Proposition 3.2. The sequence {p>} is bounded in L°°(0, T;L N ~ 2 (M JV )) when N > 3, or 
L°°(0, T; L q (M. 2 )) for all q > 1. There exists a subsequence of {pi}, still denoted by itself, 
such that 

^'(x,t)-» p(x,t), (3.20) 

strongly in C{[0,T\; L^ oc (R N )), (3 G [1,^), as j -» 00. Fere, p G L°°(0, T; L 1 n (R N )) 
when N > 3, or p £ L°°(0, T; L q (R 2 )) for all q > 1. Moreover, p(x,i) = p(r,i) is a spherically 
symmetric function. 

Proof. We only consider the case N > 3, since the proof of the case that N = 2 is similar. 
It follows from §3M) and (l3TT5|) that { yv^} is bounded in L°°(0, T; L^I^)) for 9 G [2, 2Ar 1 



JV— 2J- 

Thus, V} is bounded in L°°(0,T; (R N )) and {p?L™} is bounded in L°°(0,T;L^i 



due to (|3.12p - (|3.13p . Then, the continuity equation yields that {dtp 3 } <i R > n is bounded 

in L°°(0, T; VF -1 '^- 1 (Bi )), for any > n 2Ar . Moreover, since Vp 7 = 2^/p^V yfpi ', we have 

fe N 

{Vp } } £ < 1 #> n is bounded in L°°(0, T; L^- 1 (J3i n )). From the Aubin-Lions lemma, we get 



JV 

k 

N 

k 

Since 



^'(x,t) -> p(x,t), strongly in C([0, T]; L~ (Bi )), as j -» 00. (3.21) 



C 

Hp 3 — oil jv < -rllp J — pll + Wp 3 — p\\ jv , 

L°°([0,T];L*T=T(B„)) fc" " £°°([0,T];£Sr=2 (B i )) £°°([0,T];£W=T(Bi )) 

7« 7c'" 

we get 

^'(x,t)-»p(x,t), strongly in C([0,T]; L*=i (B n )), as j 00. (3.22) 
Clearly, ()3.20p holds and p(x, t) is spherically symmetric. □ 



10 



From (|3.9p , Lemma 13.21 and Proposition 13.21 we immediately obtain the following lemma. 
Lemma 3.3. There exists a subsequence of {p 3 }, still denoted by itself, such that 

(p0 7 (M)->P 7 (M), (3-23) 
strongly in Lj oc (R N x [0,T]), as j — > oo. 

Proposition 3.3. For any k > n 2N , there exists a subsequence of {p 3 } e <i R >n! still denoted 
by itself, such that 

vJp 3 (x,t)^V^p(x,t), W>* mL°°([0,T],L 2 (5i J), 

V k 3 

Vh(pi(x,t)) A Vh(p{x,i)), weak-* in L°°([0, T], L 2 (i?i )), 

as j -> oo, w/iere ft satisfies h(0) = and 7t'(s) = Moreover, V^/p G L°°([0, T],L 2 (R N )) 

andVh(p) e L°°([0,T],L 2 (R N )). 

Proof. It follows from (13351) that {Vy9l .< i ».>„ is bounded in L°°(0, T; L 2 (Bi )). Thus, 
there exists a function / G L°°(0, T; L 2 (Bi n )) such that, up to a subsequence, 

vJp>(x,t)^f, weak-* mL oo ([0,T],L 2 (B ln )). (3.24) 

Combining it with f)3.20[) . one can easily obtain / = V-^/p and 

\\^Vp\\l°°{[o,t],l 2 (b 1 )) < liminf ||Vv / P^||l°°([o,t],l 2 (Bi )) < C 

with a constant C independent of k and n. Clearly, we have V-^p E L°°([0, T], L^M^)). 
Similarly, we can easily obtain the result for h. □ 

From Propositions I3.1H3.2I and Corollary 13.11 using similar arguments as that in the proof of 
Lemmas 4.4 and 4.6 in [21], we can obtain the following proposition. 

Proposition 3.4. 1) Up to a subsequence, m 3 = p J U- J converges strongly in Lj oc (R N x [0,T]) 
and L 2 (0,T;Lf oc (M N )) to some m(x,t), for all (3 G [1, jfc-). 

2) \[p~ 3 V>' 3 converges strongly in Lf oc (R N x [0, T]) to -^L (defined to be zero when m = 0). In 
particular, m(x, t) = a. e. on {p(x, t) = 0} and there exists a function U(x, t) such that 

m(x,t) = p(x, i)U(x, i). 

Proof. We only consider the case N > 3, since the proof of the case that N = 2 is similar. 

1) Since {VP 7 } is bounded in L°°(0, T; L 2 f~)L^) and {^/fPV 3 } is bounded in L°°(0,T;L 2 ), 
we have that 

{p 3 V 3 } is bounded in L°°(0, T; L 1 n L^i (R N )). (3.25) 

Since V(p J U J ) = 2^/plWV '+ \J~pi ~ sffiVW , from QZJft and (l3TTl>([3TT4"l) . we obtain that 
{V(p 7 U ;? )} <i p. >n is bounded in L 2 (0, T; L 1 (-Bi )). In particular, we get 

{(p^)} <i p >„ is bounded in L 2 (0, T; W 1,:1 (jBi )). 

J — k ' J — fc ' 
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Since {p 3 } £ . < i R > n ^ s bounded in L°°(B\ x [0,T]), from (|3.2p , we can obtain that 

3 — k 5 3 — fc ' 

{WlP)} <i „ >n is bounded in L 2 (0,T;W~ 2 '^ (Bi )). 
From the Aubin-Lions lemma, we have 

p^lP — > m, 

strongly in L 2 ([0, T];L^(Bi J) for all /3 € [1, -#r). From ([535]) . we can easily obtain that 

k ' 

p^lP — ► m, 

strongly in L 2 ([0,T];L^B n )) for all /3 e [1, jfa). 

2) Using the similar argument as that in the proof of Lemma 4.6 in [21], we can obtain the 
part 2) of Proposition 13.41 where 

u = { * if ^ ' 

\ 0, ifp = 0, 

and omit the detail. □ 

Then, using similar arguments as that in the proof of Corollary 4.2 in [10], we can obtain 
the following corollary and omit the details. 

Corollary 3.2. Let mP (r, i) = {pPv?)(r,t), then 

1) there exists a function m(r,t) such that m(x, t) = m(r, t)p and m? (r,t) converges tom(r,t) 
strongly m L 2 (0, T; Lf oc ([0, oo); r^dr)) for all (3 £ [1, ^ 

2) there exits a function u(r,t) such that U(x, t) = u(r,t)~ and \fpiv? converges to ^= 
(defined to be zero when m = 0) strongly in L 2 (0, T; L 2 OC ([0, oo); r N ~ l dr)). 

Now, we show that (p, U) obtained in Propositions I3.lll3.4l satisfies the weak form of (|2.ip . 
that is (gXP holds. 

Proposition 3.5. Let (p, U) 6e t/ie Zzmii described as in Propositions \3. Then 112.11]) 

holds. Moreover, pe C([0, oo); L 1 ^)). 

Proof. We only consider the case t\ > 0, since the proof of the case that t\ = is similar. 

At first, we derive the weak form of (|3.4p . For any (p £ C*([0, oo) x [0, oo)), there exists 
n > such that supp^(-,t) C [0,n]. It follows from ([331), and (^9]) - ([3TT0|) that 

p> w N - l dr\ l l - / / + piv?<pr)r N - 1 drdt = 0, (3.26) 

Jt\ JO 

for any j satisfying Rj > n. From Proposition 13.21 we have 

OO POO 

pp ipr N ~ 1 dr — > / pipr N ~ 1 dr, 
o Jo 

and 

/•t2 /*oo />i2 /"oo 

/ / p>ip t r N - x drdt^ \ I piptr^drdt 
Jh Jo Jti Jo 
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as j — ► 0. From Proposition 13.21 and Corollary 13.21 we have 

t-2 poo pt% poo 

/ p i u j y? r r N - 1 drdt= / / \^u j )^ r r N ^ drdt -► 

t\ Jo Jti Jo 

j-t2 poo pt2 /•oo 

/ / yfp{^J~pu)Lp r r N - 1 drdt= / / puip r r N ~ l drdt, (3.27) 
Jti Jo Jti Jo 

as j — > 0. 

Therefore, taking limit j — > oo in (|3.26jl . we get 

oo pt-2 roo 

pipr*' 1 dr\ll ~ / 0"Pt + puip r )r N ' l drdt = 0. (3.28) 

Jti Jo 

For any fa G C^M^ x [ii,t 2 ]), define 

<p(r,i) = / (f>i(ry,t)dS y , 
Js 

where the integral is over the unit sphere 5 = S^ 1 in R N . Then is follows from ([3T28]) that 
([2TTT1) holds. 

Similarly, we can easily obtain 

d t ^p + div(^pu) - ~Q = 0, in V, 

where Q £ L 2 (0, T; L 2 (R N )) is the weak limit of {y^dW} in L 2 (0, T; L 2 (M N )). Thus, we 
have yfp t G L 2 ([0, oo); // _1 (lR Ar )). Since ^/p G L°°([0, oo); ff 1 ^)), we can easily get that 
VpGC([0,oo);L 2 (R iV )) . □ 

In the following, we prove that (p, U) satisfies (12.121) 

Proposition 3.6. Let (p, U) be the limit described as in Propositions \3.1]\3.4\ Then $2.12\) 
holds. 

Proof. For any (j) G C 2 ([0, oo) x [0,T]) with 0(0, t) = 4>(r,T) = 0, there exists n > such that 
supp0(-,t) C [0,n]. It follows from (f53 ]> - (|3^ and (|3T9]) - (f37T0]) that 

oo 



p u 4>(r, 0)r dr 



+ f T Hip'u^t + P>(u j ) 2 <t>r + (p*) 7 (tfr + (iV l)<t) ))r N - l drdt 
Jo Jo r 

' ^, ( (*-i)t^ + (*-i)u^ + (^-D(^-^ )tJy - ldr(ft + e i 

7 £;j r r r 2 

+ r(g(p j )+8e j ( ( S) e )(ui + {N ~ 2 1)MJ )(^ + (iV ~ 1)0 )r jV - 1 dr(it, (3.29) 



/0 ^£ 

for any j satisfying £j < i and Rj > n, where 



-b 



! I \l(2h( f S)+g(p>) + ee J (p>) e )ui}(e j ,t)ef- 1 ^ 



(3.30) 
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Claim: 



Since 



lim ^ = 0. (3.31) 



/ [(^y<t>](ej,t)dt 



T 



o 



i6[0,T] Jo 



r 



< C max ]</)(£ j, t)\ / [(^r + |d r (p>) 7 |K'- 1 drcft > 

*6[0,T] Jo ./£, 

{fpyr^drdt < C, 




Je, 




I^V') 7 !^"" 1 ^ < C / / [{fPy + \d r (p>)i\ 2 ]r N - 1 drdt < C, 

'0 Jsj Jo Jej 

and lim £j _ > o+ m a x te[o,T] l<K e j)£)l = 0, we have 

lim ef" 1 r^r^sj^dt = 0. (3.32) 

From (|3.4p and u(ej,t) = 0, we get 

(e,- ,t) + pP (sj , t)d r u j (ej , t) = 0. 

Thus, using (|2.4jl . we have 



T 



2hUP)+g{pP) 



ef" 1 jf {(2h( fi >)+g(jP))vi<l>){e j ,t)dt = -e* r - 1 J f ' ' ' ^^o)( -,,./)>// 
= -ef- 1 J\2d t h(pj)0)(e v t)dt 

= 2ef- 1 hU(e j ))cf>(e j ,0) + fcf" 1 jf* \h(pP)d t 4>){e j ,t)dt. 
It is easy to obtain 



JV-l 
2 

'J 



\y/p^t)\ < C\\^\\ HH[E . iRj]) < Ce] 

\\[7^)\<c\\^\\ miXxm <c, 

JV-l 

< C|M^'(l,t))| + ||V^|| i2([£f)1D < C + C £j 2 , (3.33) 
where /i satisfies /i(0) = and hf(s) = Since h(s) < yfsh^s), we have 

Thus, we can easily obtain 



lim max \dt<b(Ei,€)\ =0 



and 



(h(p>)dt<j>)(ej,t)dt 



< C max |9{^(e 7 -,t)| — > 0, 
te[o,T] 
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as Ej — > 0+. Hence, we have 



rT 

JV-1 



lim e 1 ]-' \ {{2h{p>) + g{p>))ul<t>){e j ,t)dt = Q. 



o 



Similarly, one can obtain that 

lim sf [ T [(pi) 9 44>}(e 3 ,t)dt 
£ j^ u + Jo 

Thus, dSM]) holds. 



Now, for any fa G (C%(R N x [0,3^))^ with supp0 2 (-,t) C B n and </> 2 (x,T) = 0, we set 

cf>(r,t) = ( Mry,t)-ydS y . (3.34) 



s 

Since 

(r JV - 1 0)r = S r / div0 2 (x,t)(ix = / (^.(n/,*)^, 

J|x|<r 

we have by direct calculation that 

rT poo (at -i\„.iA rT 



-[ 2h(p>)(uUr + (iV l^^ y^drdt = - [ [ 2h{p>)D{W):V<p 2 dxdt. 

JO Jej r JO J\x\>Ej 

Similarly, one has 



J Sj r r 

I W) + eej(fjf)<iivW&w<t>2dx.dt 

J\x\>Ej 



and 



[ (p>) e (divU j div<j>2 - D(U j ) : Vfo)dxdt. 

J\x\>ej 

Thus, from f)3.29p . we have 

f /9 J U J • 2 (x, 0)dx + f [ (\/^y^\J j ■ d t (f>2 + y^\J j : Vfo)dxdt 

J\x\>ej Jo Jm. n 

rp rp 

(piydivfcdxdt- f [ [2h{p>)D(XJ j ) : Vfa + g(p>)divV j div fa]dxdt 
o Jr n Jo J\x\>ej 



Ej ft [(6- l)(^ydivU J div0 2 + (p>) 9 D{W) : Vfojdxdt + e{. (3.35) 

JO J\x\>ei 



We proceed to show that each term on the left hand side of (|3.35p converges to corresponding 
term in (|2.12p . and each term on the right hand side of f)3.35j) vanishes as j — ► oo. 
First, the proof of the convergence of p'XJWtfa is similar to that of (13. 27ft . 
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Next, from Proposition 13.41 we obtain 
rT 



piW ® VpJ'U J : V<p 2 dxdt 



/0 JR N 
From Lemma 13.31 we have 



o Jr n 



y^U (g) y'joU : V(p 2 dxdt, as j — > oo. 



JM^ 



(/7 7 ) 7 div02C?xdt 



JR N 



p 1 div (f) 2 dxdt, as j — > oo. 



Concerning the diffusion terms on the left hand side of (I3.35p . using (|3.9I) and integration by 
parts, we have 



2h{p>)D{W) : Vfodxdt 



J\x\>£j 

T 



JR N 
T 



^l(^W) ■ Afo + ^H^W) • Vdiv0 2 
VP 3 v P 3 



dxdt 



[(VP j V J ) ■ (Vh(p>) ■ V)0 2 + (y/pHJ 3 ) • (Vfc ■ V)h(pP)]dxdt. (3.36) 



'0 JB e ., n 

Using the similar argument as that in the proof of (|3,33p , we have 



< C n , and Hp^llioo^i^jxjo^]) < Cfe,n- 

L°°([0,T];L 2 (B n )) 

Then, using the similar argument as that in the proof of (|3.27p . we have 



and 



Hp 3 ) , 

./r^ v//^ 



^IP) • kfadxdt 



f (vV'U 5 ) • Vdiv^dxcft 



Jr n \[P 

Hp) 



(y/pU) ■ Afodxdt, (3.37) 



o Jr n \ff> 



lo JR N y p> 

as j — > oo. From Corollary 13.11 Lemma 13.21 and Propositions I3.3 l l3.4l we have 
rT 



(VJoU) • Vdiv^dxcft, (3.38) 





T 



V'lP) • (yh(p>)-V)<fadxdt 



< 



o is in{|UJ'|<Af} 



/ 

JB i n 



{|UJ|>M} 



( WlF) • (yh(p>) ■ V)fodxdt 



7(iV+l)-jV 

v ' ' ; \ L W (R^x[0,T]) * 



+ 



1 



1 + ln(l + M 2 ) 

0, as M, — » oo, 



|HU^(l + ln(l + |U^))|| Loo(0iT;Ll) 



r 







(VpU) • (V/i(p) • V)cp 2 dxdt ->■ 0, as fc -> oo, 
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and 



piW) ■ (Vh(p>) ■ V)<fadxdt 



JBx 

k ' 



(y/pU) ■ (V/i(p) • V)4> 2 dxdt 



JBx 

k ' 



as j —* oo. Thus, we have 



WpiW) ■ {Vh{p>) ■ V)4> 2 dxdt 



JB ej , n 



as j — > oo. Similarly, we can obtain 







(VpU) • (V/i(p) • V)4> 2 dxdt, (3.39) 



(V>'U J ) • (V0 2 • V)h(p>)dxdt 



J5 Ej ,„ 



as j — » oo. From (|3.36p - (|3.40p . we obtain 

'|x|> £j 

Similarly, we obtain 



o ^i" 



(vW • ( v< ^2 • V)h(p)dxdt, (3.40) 



/ / 2/i(/>?)£>(U j ) : Vfodxdt -►< 2h(p)D(U), V<£ 2 >, as j -» oo. 



g(p J )divU J div(/)2(ixdt — >< g(p)divU, div</> 2 >, as j — > oo. 

^|x|>£j 

Up to now, we have proved that each term on the left hand side of (|3.35p converges to 
corresponding term in (|2.12p as j — > oo. In the following, we prove that each term on the right 
hand side of (|3.35p vanishes as j — > oo. 

From Lemma |3.1[ we get 



JR N 



(piydivWdivfcdxdt 



< cyeJll v &IUs( £ . 

Nii-e) 

< C{T)^e]n—^- 



(p>Y\VV J \ 2 dxdt) 

JR N J 



T \ I 

/ / p j dxdt) \B n \^- 
Jo Jr n ) 



and 



N(i-e) 

<C(T)^e]n—. 



(3.41) 
(3.42) 



(p>)°D(U 3 ) : Vfodxdt 

It follows from (|3.3ip and (|3.4ip - (|3.42p that each term on the right hand side of (|3.35p vanishes 
as j — > oo. 

Taking the limit j — > oo in (|3.35|) . we finish the proof of this proposition. □ 



From the above arguments, we can immediately finish the proof of Theorem 12.11 
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